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On a Density for Sets of Integers 
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Abstract 

A relationship between the Riemann zeta function and a density on integer sets is explored. Several properties of the examined 
density are derived. 
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1 Introduction 


of A is given by 


Several measures for the density of sets have been discussed in 
the literature ma. Presumably the most employed tool for eval¬ 
uating the density of sets is the asymptotic density, also referred 
to as natural density. The asymptotic density is expressed by 


d(A) = lim 

n—>°o 


l|An{l,2,3,...,»}|| 

n 


(1) 


£d(A) = lim ^ a '- n a - . (5) 

n^co log/? 

InlH, Bell and Burris bring a good exposition on the Dirich- 
let density. The Dirichlet density is defined as the limit of the 
ratio between two Dirichlet series. Let A C B be two sets. The 
generating series of A is given by 


provided that such a limit does exist. The symbol || • || denotes 
the cardinality, and A is a set of integers. Analogously, the lower 
and upper asymptotic densities are defined by 


A(s) 


E 

n— 1 


N(A,n) 


( 6 ) 


d(A) = liminf 

n 


d(A) = lim sup 


l|A l~l {1,2,3,...,«} 


( 2 ) 

(3) 


The asymptotic density is said to exist if and only if both the 
lower and upper asymptotic densities do exist and are equal. 

Although the asymptotic density does not always exist, 
the Schnirelmann density |2][3j is always well-defined. The 
Schnirelmann density is defined as 


5(A) = inf 

n> 1 


]|A n {1,2,3,... ,«}H 
n 


(4) 


Interestingly, this density is highly sensitive to the initial ele¬ 
ments of sequence. For instance, if 1 ^ A, then 5(A) = 0 |7j. 

Another interesting tool is the logarithmic density fflO- Let 
A = { 01 , 02 , 03 ,...} be a set of integers. The logarithmic density 
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where N(A,n) is a counting function that returns the number of 
elements in A of norm n (T). The Dirichlet density is then ex¬ 
pressed by 


D(A) = lim 


A(s) 


sfa Bfs) 


(7) 


where B (s) is the generating series of B and a is an abscissa of 
convergence Q. In l6j], we also find the Dirichlet density defined 
as 

D (A) =lim(s- 1) £ 2 (8) 

whenever the limit exists. This density admits lower and up¬ 
per versions, simply by replacing the above limit by lim inf and 
lim sup, respectively. 

The aim of this work is to investigate the properties of the 
Dirichlet density as defined in Equation [7] in the particular case 
where the set B is the set of natural numbers. This induces a 
density based on the Riemann zeta function. 


2 A Density for Sets of Integers 

In this section, we investigate the particular case of the Dirich¬ 
let function, applied for subsets of the natural numbers. In this 
case, taking into consideration the usual norm, where the norm 
of any natural number is equal to its absolute value, the counting 
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function becomes 


Proof: dens(A c ) = dens(N— A) = dens(N) — dens(NHA) = 1 — 


. f 1, if n G A, 

N(A,n) = < ’ 

10, otherwise. 

(9) 

dens (A). 

Proposition 2 (Monotonicity) The discussed density, 
tone function, i.e., if A C B, then dens(A) < dens(B). 

is a mono- 

Definition 1 The density dens of a subset A C N is 

given by 

Proof: We have that 


dens(A)4lim I ' r=lA?(A, " ) Y 

4i E“ =1 lV(N,n)£ 

= lim L,6A t 

41 y°° ± 

v L*n= 1 n s 

(10) 

dens(B) = lim 

41 CM 

(19) 

(11) 

= lim ^ ,,eA F + lim ^ neB ~ A F 

41 CM 41 CM 

(20) 

1 • ^*neA fjs 1 

= hm , j > 1, 

41 Co) 

(12) 

li EneAF 

-41 CM 

(21) 

if the limit exists. The quantity £(•) denotes the Riemann zeta 
function fffj. 

= dens (A). 

(22) 

□ 

Proposition 1 The following assertions hold true: 


Proposition 3 (Finite Sets) Every finite subset of N has density 

1. dens(N) = 1 


zero. 


2. dens(A) > 0 (nonnegativity) 


Proof: Let A be a finite set. For .v > 1, it yields 


3. i/dens(A) and dens(B) exis^and A flB = 0, then dens(A U 
B) exists and is equal to 

> — < > - = constant < 00. 

^ yiS ^ yi 

neA n neA n 

(23) 

dens(A) +dens(B) (additivity). 


Thus, 


Proof: Statements 1 and 2 cab be trivially checked. The additiv¬ 
ity property can be derived as follows: 

<i« s (A) = li m Z: ”fy<Um CO “ 1 “ , =o. 

41 CM 41 CM 

(24) 

dens (A UB) = lim »' 

41 CM 

(13) 

□ 

As a consequence, sets A C N of nonzero density must be infinite. 

= lim ^ neA W + lim eB 1 

41 CM 41 CM 

= dens(A) +dens(B). 

I 

II (14) 

Corollary 4 The density of a singleton is zero. 


(15) 

Proposition 4 (Union) Let A and B be two sets of integers. Then 
the density of A UB is given by 


□ 

dens(A UB) = dens(A) +dens(B) — dens(A FlB) 

(25) 


Corollary 1 The density of the null set is zero. 

Proof: In fact, dens(0) =0, since 1 =dens(N) =dens(0UN) = 
dens(0) +dens(N) = dens(0) + 1. □ 

Corollary 2 dens(B — A) = dens(B) — dens(A flB). 


Proof: The proof is straightforward: 

dens(B) = dens((A ITB) U (A e flB)) (16) 

= dens(AnB)+dens(A c nB) (17) 

= dens(A flB) + dens(B — A), (18) 

since A (IB and A c (IB are disjoint. □ 


Corollary 3 For every A, dcns(A') = 1 — dens(A), where A c is 
the complement of A. 

1 l or ease of exposition, in the following results, we assume that the densities 
of the relevant sets always do exist. 


Proof: Observe that AUB = AU(B- A), and A f] (B — A) = 0. 
Then, it follows directly from the properties of the density mea¬ 
sure that 

dens(A LIB) = dens(A U (B — A)) (26) 

= dens(A)+dens(B — A) (27) 

= dens(A) +dens(B) — dens(A flB). (28) 

□ 

Proposition 5 (Heavy Tail) Let A = 

{ai,ci 2 , ■ ■ ■ ,aN-haN,aN+i, ■ ■ •}• If A\ = {a\,ci 2 , ■ ■ ■ ,aiv-i} 
and Ai = {qn,gn+ i, ...}, then 

dens(A) = dens(A2). (29) 

Proof: Observe that A = A j U A9 and A] and A2 are disjoint. 
Therefore, dens(A) = dens(Ai) + dens(A2). Since Ai is a finite 
set, dens (A 1) = 0. □ 
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Now consider the following operation m <g) A = { ma \ a £ 
A,m £ N}. This can be interpreted as a dilation operation on 
the elements of A. 

In Erdos et al. examined the density of the set of 

multiples m®A , showing the existence of a logarithmic density 
equal to its lower asymptotic density. Herein we investigate fur¬ 
ther this matter, evaluating the density of sets of multiples. 

Proposition 6 (Dilation) Let A be a set, such as dens(A) > 0. 
Then 

dens(m(g>A) = — dens(A). (30) 

m 

Proof: This result follows directly from the definition of the dis¬ 
cussed density: 


Taking the limit as s j, 1, the above upper bound becomes 


lim 2 n = dens(A^) = dens (A). 

41 CO) 


( 41 ) 


Now examining the lower bound and using the dilation property, 
we have that 


lim 

41 


-neA% ' 


«N+ l . 


CO 


UN 

ClN H - 1 
GN 

ClN H - 1 


dens(A2 ? ) 
dens (A). 


(42) 

(43) 

(44) 


densfm 0 A) = lim-—-— (31) 

41 CO) 

= l im ^I^AF (32) 

41 CO) 

= — dens (A). (33) 

m 

□ 

Let A © m = {a+m \ a £ A,m £ N}. This process is called a 
translation of A by m units flOl p.49]. Our aim is to show that 
the discussed density is translation invariant, i.e., dens(A ffim) = 
dens (A), m > 0. Before that we need the following lemma. 

Lemma 1 (Unitary Translation) Let A be a set, such as 
dens (A) > 0. Then 


Since Ve > 0, 3N such that 

—dens(A) > dens(A) — e, 
aN+ 1 


and as 

lim ”4, ( ? +1 ^ = dens(Af ® 1) = dens(A® 1), 

41 CO) 

it follows that Ve > 0 we have that 

dens(A) e < dens(Affi 1) < dens(A). 
Therefore, letting e —>• 0, it follows that 
dens(A ffi 1) = dens(A). 


(45) 


(46) 


(47) 


(48) 


dens(A ffi 1) = dens(A). (34) 


□ 


Proof: Let A = {ai,a 2 , ■ ■ ■ ,aN-i,aw,aN+i ,•••}■ We can split A 
into two disjoint sets as shown below: 


Proposition 7 (Translation Invariance) Let A be a set, such as 
dens (A) > 0. Then 


A = {a\,a2, ■ ■ ■ ,a N -i,aN,a N+ i,...} (35) 

= {ai,a2,--- 1 a N -\}U{aN,a N+ \,.. (36) 
= Ai UA 2 , (37) 

where A^ is a finite set and A C is a ‘tail’ set starting at the element 
on- By the Heavy Tail property, we have that VA f dens (A) = 
dens (A 2 ). For the same reason, V/V, dens(Affil) = dens (A^ ffi 1). 

Note that, for s > 1, if kn > n + 1, then the following inequal¬ 
ities hold: 


^ neA 2 ( kn) s ^-‘neA% (n+1) 1 ^ ^ neA% n s 

CO) - CW -^0T 


(38) 


In other words, for the above inequality to be valid, we must have 


k> 


n + 1 
n 


Vrc £ A%. 


(39) 


For instance, let k = (ax + l)/«v. Consequently, we establish 
that 



CO) 


- CO) 


^ neAy n s 

-“cor 


(40) 


dens(A ffi/n) = dens(A), (49) 

where m is a positive integer. 

Proof: The proof follows by finite induction. We have already 
proven that dens(A ffi 1) = dens(A). Therefore, we have that 

dens(Affi (m+ 1)) = dens((A ffi/n) ffi 1) (50) 

= dens(Affim) = dens(A). (51) 

□ 

One possible application for a density on a set of natural num¬ 
bers A is to interpret it as the chance of choosing a natural num¬ 
ber in A when all natural numbers are equally likely to be chosen. 
Interestingly, as we show next the above measure of uncertainty 
does not obey all the axioms of Kolmogorov since it is not a- 
additive. Additionally, we show that it is impossible to define a 
finite cr-additive translation invariant measure on (N,2 n ). This 
result emphasizes an important point that there are reasonable 
measures of uncertainty that do not satisfy the formal standard 
definition of a probability measure. 
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Theorem 1 There is no a-additive measure defined on the mea¬ 
surable space (N,2 n ) such that: 























1. 0 < jU(N) < oo; and 

2. p is translation invariant. 

Proof: Suppose that p is translation invariant. Then every single- 
ton set must have the same measure. Let ©i < fth be any natural 
numbers. Since {C 0 i} ® (ft) 2 — © 1 ) = {© 2 }. we have 

M({®2}) = fi({®i}©(©2-©i)) =At({©t}), 

where the last equality follows from translation invariance. Let 
C = ju({1}). If c = 0, then p{N) ^0 = Z?=iP ({«'})• Thus p is 
not cr-additive. If c > 0, then fi(N) < °° = p ({;}), which 
also implies that p is not <7-additive. □ 

Proposition 8 (Criterion for Zero Density) Let A = 

{a l ,a 2 ,...,a„,a n+1 ,...}. If lim 41 £r=i 3 j converges, then 
dens (A) = 0. 


Corollary 7 (Geometric Progressions) Let G = 

{ar,ar 2 ,ar 3 ,...}, where a and r > 1 are positive integers. 
Then 

dens(G) = 0. (55) 

Proof: Notice that, for r > 1, 

= I£J-<I£ (T)"<oo (56) 

„=i ( ar "Y aS „=i rm a$ n =i W 

Thus, it also follows from the preceding discussion that 

dens(G)=0. □ 

Let M p = {p,2p,3p ,...} be an arithmetic progression, where 
p G N. For example, = {2,4,6,...}, M 5 = {5,10,15,...}, etc. 
Regarding the cardinality of these sets, we have that \\M p \\ = Kg. 

Proposition 9 (Arithmetic Progressions) For a fixed integer p, 
the density of the set M p is given by 


Proof: It follows directly from the definition of dens (A). □ 

Definition 2 (Sparse Set) A zero density set is said to be a 
sparse set. 

As matter of fact, any criterion that ensure the convergence of 
lim 4 i fff | /- can be taken into consideration. In the next result, 
we utilize the Ratio Test for convergence of series. 


dens(Mp) = 


(57) 


Proof: Note that M p = p Cg LI. Thus, according to Proposition [ 6 ] 
we have that 


dens(Mp) = dens (pig) N) = — dens(N) = —. 


(58) 


Corollary 5 Let A = {cq ,a 2 , ■ ■ ■ ,a Hl a „ + \,...}. If 

lim^oo < 1, then dens (A) = 0. 

a n+1 

Proof: Observe that for s > 1 



According to the Ratio Test El p. 68 ], the series on the right- 
hand side of the above inequality converges whenever 

1 

lim = ii m —— < 1 ( 53 ) 

«->°° — n-s-oo a n +\ 


□ 

Consequently, we have, for instance, dens(Mi) = dens(N) = 1, 
dens(M 2 ) = 1/2, dens(M 5 ) = 1/5 etc. 

One can derive a physical interpretation for the density of 
M p . Consider discrete-time signals M p [n] characterized by a se¬ 
quence of discrete-time impulses associated to the sets M p . The 
signals M p [n] are built according to a binary function that returns 
one if n is a multiple of p, and zero otherwise. For example, we 
can have M 2 [n\ and M^[n\, associated to M 3 and M 5 , as depicted 
in FigureQ] In terms of signal analysis, dens(M p ) has some cor¬ 
respondence to the average value of M p [n\. 

3 Density of Particular Sets 


Thus, by series dominance, it follows that Y.neA "v also con¬ 
verges. And finally, this implies that the density of A must be 
zero. □ 


In this section, some special sets have their density examined. 
We focus our attention in three notable sets: (i) set of prime num¬ 
bers, (ii) Fibonacci sequence, and (iii) set of square-free integers. 


Lemma 2 (Powers of the Set Elements) Let A = {a u a 2 ,...} 3.1 Set of Prime Numbers 

and m > 1 be an integer. The density of A m = {a'! 1 ,a'f,...} is 

zero " Let P = {2,3,5,7,11,13,17,19,23,...} be the set of prime num¬ 

bers. 


Proof: Notice that, for ,v > 1 and m > 1, 


tt «) 


1 f 1 f 1 

^ L nm ^ ^ ' 


(54) 


Thus, it follows from the preceding discussion that dens(A) = 0. 

□ 


Corollary 6 The density of the set of perfect squares is zero. 


Proposition 10 The prime numbers set is sparse. 


Proof: Utilizing the definition of the prime zeta function, C,'(s) = 
Lr.pi is primep/ 5 ’ we have that den s (P) = hniq, Addition¬ 

ally observe that f,'(s) < °° for s > 1. Taking in account that 


In COO 



(59) 
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1.2 



n 


(a) 



n 

(b) 

Figure 1: Discrete-time signals M 3 [n] and M 5 [n] associated to the 
sets M 3 and M 5 , respectively. 


1.2 



n 

Figure 2: Discrete-time signal P[n] corresponding to the set of 
prime numbers. 

Since the density is zero, it indicates that the associated 
discrete-time signal P[n\, as shown in Figure [3 has null average 
value. 

3.2 Fibonacci Sequence 

The Fibonacci sequence constitutes another interesting subject 
of investigation. Fibonacci numbers are constructed according 
to the following recursive equation u[n] = u[n — 1] + u[n — 2], 
n > 1, for m[0] = u[ 1] = 1, where u[n\ denotes the nth Fibonacci 
number lfl2l p. 1 60] . This procedure results in the Fibonacci set 
F = {1,2,3,5,8, 13,21,34,55,.. .}. 

Proposition 11 The Fibonacci set is sparse. 


for e > 0 we have that 


lnC(l+e) = E 


cm+e)) 


k= 1 


= ?'(l + e) + I 

k—2 

Now consider the following inequalities: 


C'(*(l + e)) 


0 < 


C'(l+e)_lnC(l + e)-ir=2 m ^ ±£l1 


C(l + e) 


< 


C(l + e) 

lnC(l + e) y M±£ll 

C(l+e) hw+ ^ 

lnC(l+£) 


C(l + e) ‘ 

As e -A 0, we have that C,{\ + e) —>• 00 , therefore 

lim ln j (1 +£) = 0, 


(60) 

( 61 ) 

(62) 

(63) 

(64) 

(65) 


C(l + e) 

because lim x _>ooln(^)/a: = 0. Therefore, we can apply the 
Squeeze Theorem once more, and find that 


C (1 +e) 


( 66 ) 

□ 


Proof: It is known that the sum of the reciprocals of the Fi¬ 
bonacci numbers converges to a constant (reciprocal Fibonacci 
constant), whose value is approximately 3.35988566... El. 
Consequently, for s > 1, we have that 

^Ldn —2 u\ti\ s 

dens(F) = lim— 

41 C(«) 

I“=i i 3.35988566 

<llm —? 7 T- L = h m -77T— 

41 C(s) 41 C(^) 

□ 


(67) 

— = 0 . ( 68 ) 


3.3 Set of Square-free Integers 

The set of square-free numbers can be defined as S = {n £ 
N | |jLt(«)| = 1}, where ju.(•) is the Mobius function, which is 
given by 


p(n) 


1 , if n = 1 . 

0 , if p 2 \n, for some prime number p. 

(—l) r , if n is the product of distinct 
prime numbers. 


(69) 


Thus the first elements of S are given by 
1,2,3,5,6,7,10,11,13,14,15,... 
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Proposition 12 The density of the set of square-free integers is 

WY 

Proof: It is known fT4l that £“ =1 thus it follows 

easily that: 


dens (S) = lim 

41 C(s) 

(70) 

y°° ImMI 

= lim - " 

41 £(s) 

(71) 

CM 

= lim= lim ■... , 

41 C(s) 41 C(2s) 

(72) 

1 6 
_ C(2) _ it 2 ' 

(73) 


□ 



Figure 3: Computational result for the approximate density 
for the truncated prime set Mf, according to (i) the function 
dens' (P T ) (bold curve), and (ii) the function d '(P T ) (thin line). 
Observe the convergence to 1/2. 


4 Computational Results 


Definition 3 Let A be a set of positive integers. The truncated 
set AJ is defined by 

A r =An{i,2,3,..,r-i,r}, (74) 


where T is a positive integer. 

In other words, A r contains the elements of A which are no 
greater than T. 

The concept of truncated sets can furnish computational ap¬ 
proximations for the numerical value of some densities. Con¬ 
sider, for example, the quantity 


d'(A r ) 


A Ml 

iiN r ir 


(75) 


This expression d'(-) can be interpreted in a frequentist way as 
the ratio between favorable cases and possible cases. Clearly, the 
asymptotic density can be expressed in terms of d'(-): 


Definition 4 (Approximate Density) The approximate density 
for a truncated set AJ is given by 

dens \A r ) = ^^f. (80) 

Lm=1 n 

Again, whenever the limit order of Equation [77] can be inter¬ 
changed, we have that 

dens(A) = lim dens'(A r ). (81) 

T—> °o 

The quantity dens'(-) furnishes a computationally feasible way 
to investigate the behavior of dens(-). 

In the following, we obtain computational approximations for 
the asymptotic density and the discussed density. The approx¬ 
imate densities are numerically evaluated as T increases in the 
range from 1 to 1000. Now we investigate (i) arithmetic pro¬ 
gressions, (ii) the set of prime numbers, and (iii) the Fibonacci 
sequence. 

4.1 Arithmetic Progressions 


d(A) = lim d'(A T ). (76) 

In a similar fashion, we can consider a version of the discussed 
density for truncated sets as defined below. First, observe that 
discussed density dens(-) can be expressed as the following dou¬ 
ble limit 


V — 

dens (A) = lim lim . (77) 

v ’ 41 r->°o Yj L 

Y L^n = 1 n s 

Restricted to the class of sets for which the above limits can have 
their order interchanged, we find that 


dens (A) = lim lim 

T ^oo 


£ neA T 

~yf T~ 

Zw/7 — 1 tl S 


( E«eA r j; 

V iLl ~n 


(78) 


Considering an arithmetic progressions M q , we have that: 

d(M q )= limd'(M[) = l (82) 

T—f oo ^ q 

and 

dens (Mq )= lim dens' (Mj). (83) 

Both quantities limits converge to the same quantity 1 /q. Intu¬ 
itively, we have that the chance of “selecting” a multiple of q 
among all integers is 1 /q. This is exactly the density of M q . Fig¬ 
ure [3] shows the result of computational calculations of d'(-) and 
dens'(•) for Mf. 

4.2 Prime numbers sets 


= lim 

T—>oo 


(79) 
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Consider the truncated set of prime numbers P T = {p : 
p is a prime,/? < T). One may easily compute both d \P T ) and 
dens' (P T ) for any finite T. 




















Figure 4: Computational result for the approximate density 
for the truncated prime set P T , according to (i) the function 
dens' (P T ) (bold curve), (ii) the function d'(P T ) (thin line), and 
(iii) the approximation based on the Prime Number Theorem 
(dotted line). 


An alternative path for the estimation of d '(P T ) = where 

the function n(x) represents the number of primes that do not 
exceed x, is to utilize the approximation for K ( T) given by the 
Prime Number Theorem 11211151 p.336]: 

*(r)«ii(r), (84) 

where li(-) denotes the logarithmic integral lfl6l . The curves dis¬ 
played in Figure [4] correspond to the calculation of dens'(P 7 ), 
d '{P T ) and li(7') jT. As expected, all curves decay to zero. 


4.3 Fibonacci Sequence 

Now consider the truncated Fibonacci set F T = 
{1,2,3,5,8,13,21,... ,r}. The result of calculating 


and 


d'(P r ) 


El 

l|N r || 


dens' (F t ) = 


^ neF T n 

~yT F 

J-n=l n 


(85) 


( 86 ) 


are shown in Figure[5] Both curves tend to zero as T grows. This 
fact is expected since we have already verified that dens(F) = 0. 
The small convergence rate of dens'(•) is typical of problems that 
deal with the harmonic series. 


5 Conclusion 

A density for infinite sets of integers was investigated. It was 
shown that the suggested density shares several properties of 
usual probability. A computational simulation of the proposed 
approximate density for finite sets was addressed. An open topic 
is the characterization of sets for which the limiting value of the 
approximate density equals their density. 


Figure 5: Density of the truncated Fibonacci set: dens'(P r ) 

(bold line) and d \F t ) (thin line). 
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